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We present a detailed analysis of the zenith angle distributions of atmospheric neutrino events 
observed in the Super-Kamiokande (SK) underground experiment, assuming two-flavor and three- 
flavor oscillations (with one dominant mass scale) among active neutrinos. In particular, we calculate 
the five angular distributions associated to sub-GeV and multi-GeV /r-like and e-like events and to 
upward through-going muons, for a total of 30 accurately computed observables (zenith bins). First 
we study how such observables vary with the oscillation parameters, and then we perform a fit to the 
experimental data as measured in SK for an exposure of 33 kTy (535 days). In the two-flavor mixing 
case, we confirm the results of the SK Collaboration analysis, namely, that <-> u T oscillations are 
preferred over zy, <-> v e , and that the no oscillation case is excluded with high confidence. In 
the three-flavor mixing case, we perform our analysis with and without the additional constraints 
imposed by the CHOOZ reactor experiment. In both cases, the analysis favors a dominance of the 
i'll <-> v T channel. Without the CHOOZ constraints, the amplitudes of the subdominant v M <-► z/ e 
and v & <-> v T transitions can also be relatively large, indicating that, at present, current SK data 
do not exclude sizable v e mixing by themselves. After combining the CHOOZ and SK data, the 
amplitudes of the subdominant transitions are constrained to be smaller, but they can still play a 
nonnegligible role both in atmospheric and other neutrino oscillation searches. In particular, we 
find that the v e appearance probability expected in long baseline experiments can reach the testable 
level of ~ 15%. We also discuss earth matter effects, theoretical uncertainties, and various aspects 
of the statistical analysis. 


PACS munber(s): 14.60.Pq, 13.15.+g, 95.85.Ry 


I. INTRODUCTION 


The Super-Kamiokande (SK) water-Cherenkov experiment |jj has recently confirmed @0, with high statistical 
significance, the anomalous flavor composition of the observed atmospheric neutrino flux, as compared with theoretical 
expectations ||[|. The flavor anomaly had been previously found in Kamiokande |t]|| and IMB ||, and later in 
Soudan2 0- but not in the low-statistics experiments NUSEX |ll| and Freius 0- 

The recent SK data have also confirmed earlier Kamiokande indications Q for a dependence of the flavor anomaly 
on the lepton zenith angle which is correlated with the neutrino pathlength. The features of this dependence 
are consistent with the hypothesis of neutrino oscillations, which represents the most natural (and perhaps exclusive) 
explanation of the data Q. The oscillation hypothesis is also consistent with other recent atmospheric neutrino 
data, namely, the finalized sample of Kamiokande upward-going muons ||l3| , the latest muon and electron data from 
Soudan2 0 - and the samples of stopping and through-going muons in MACRO jl5| . 

The SK atmospheric v measurements, which are described in detail in several papers 011 , conference proceedings 
0S and theses (00, demand the greatest attention, not only for their intrinsic importance, but also for their 
interplay with other oscillation searches, including solar v experiments 123 and long baseline oscillation experiments at 
reactors 0 0 and accelerators (27-31 . In this work, we contribute to these topics by performing a comprehensive, 
quantitative, and accurate study of the SK atmospheric v data in the hypothesis of three-flavor mixing among active 
neutrinos. We also include, within the same framework, the recent data from the CHOOZ reactor experiment 10 . 0 . 

More precisely, we consider the five SK angular distributions associated to sub-GeV and multi-GeV //-like and 
e-like events j|g] and to upward through-going muons 00], for a total of 30 accurately computed observables 
(5+5+5+5+10 zenith bins). First we study how such observables vary with the oscillation parameters, and then we 
fit them to the experimental data as measured in SK for an exposure of 33 kTy (535 days) j|0]. 

In the two-flavor mixing case, we confirm the results of the SK Collaboration analysis, namely, that v T 

oscillations are preferred over v e , and that the no oscillation case is excluded with high confidence. In the three- 

flavor mixing case, we perform our analysis with and without the additional constraints imposed by the CHOOZ 
reactor experiment. In both cases, the analysis favors a dominance of the <-> v T channel. Without the CHOOZ 
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constraints, the amplitudes of the subdominant y,, 


and y e <-> v T transitions can also be relatively large, indicating 


that, at present, current SK data do not exclude sizable y e mixing by themselves. After combining the CHOOZ and 
SK data, the amplitudes of the subdominant transitions become smaller, but we show that they can still play a 
nonnegligible role both in atmospheric, solar, and long baseline laboratory experiments. 

The plan of our paper is as follows. In Section II we discuss the 30 SK observables used in the analysis, as well 
as the CHOOZ measurement. In Section III we set the notation for our three-flavor oscillation framework. The two- 
flavor subcases are studied in Sec. IV. In Sec. V we perform the three-flavor analysis of SK data (with and without 
the CHOOZ constraints) and discuss the results, especially those concerning v e mixing. In Sec. VI we study the 
implications of our analysis for the neutrino oscillation phenomenology. We conclude our work in Sec. VII, and devote 
Appendixes A and B to the discussion of technical details related to our calculations and to the statistical analysis. 

Some of our previous results on two- and three-flavor oscillations of solar [|32] Q , atmospheric |l5] [l 8 |, laboratory 
p9| , [Io| neutrino experiments and their combinations 141 - 441 will be often referred to in this work. In particular, the 
3 v analyses in ]32|,^6] [39| summarize the pre-SK situation. However, we have tried to keep this paper as self-contained 
as possible. 


II. EXPECTATIONS AND DATA 

In this Section we discuss expectations and data for five SK zenith distributions: sub-GeV e-like and /rt-like, multi- 
GeV e-like and /x-like, and upward-going muons, with emphasis on some critical aspects of both theory and data that 
are often neglected. Finally, we discuss the CHOOZ reactor results. 


A. Zenith distributions of neutrinos and leptons 

A basic ingredient of any theoretical calculation or MonteCarlo simulation of atmospheric v event rates is the flux 
d> of atmospheric neutrinos and antineutrinos 1 ] as a function of the energy E v and of the zenith angle O. The flux 
<fr(E v , O) is unobservable in itself, and what is measured is the distribution of leptons l = /i, e after y interactions, as 
a function of the lepton energy Eg (< E v ) and lepton zenith angle 9 ( 7 ^ 0). 

In Fig. 1 we show the sum of theoretical v and y fluxes, as a function of the neutrino zenith angle 0, for selected 
values of the energy E u (1, 10, and 100 GeV). The fluxes refer to the calculations of [|| (HKKM’95, solid lines) and 
0 (AGLS’96, dots) without geomagnetic corrections (so that the sign of cos© is irrelevant in this figure). The upper 
and middle panels refer to muon and electron neutrinos, respectively, while the lower panel shows their ratio. Several 
interesting things can be learned from this figure. For instance, the often quoted value + v l _ b )/(y e + V e ) — 2 for the 
muon-to-electron neutrino ratio clearly holds only for low-energy, horizontal neutrinos. This ratio increases rapidly 
as the neutrino energy increases and as its direction approaches the vertical. In fact, both y,, and v e fluxes decrease 
towards the vertical (see upper and middle panels), where the slanted depth in the atmosphere is reduced; however, 
z/ e ’s are more effectively suppressed than yfls, due to their different parent decay chains. In addition, the greater the 
energy of the parents, the longer the decay lengths, the stronger the dependence of the y^/ve flux ratio on the slanted 
depth and thus on the zenith angle 0. In other words, high-energy, vertical “atmospheric v beams” are richer in zx^’s 
and, therefore, are best suited in searches for zx M —> v e oscillations, where initial zx e ’s represent the “background” (see 
also |Q). We anticipate that, in fact, multi-GeV data are more effective than lower-energy (sub-GeV) data in placing 
bounds on —> v e transitions. Another consequence of the non-flat y^/v e ratio is the appearance of distortions of 
the zenith distributions that, although related to vacuum neutrino oscillation, do not depend on neutrino pathlength- 
to-energy ratio L/E v J38|. Finally, notice in Fig. 1 that the good agreement at low energies between the two reported 
calculations of v^/ve is somewhat spoiled at high energies. This shows that the often-quoted uncertainty of ±5% for 
the y^/ve ratio, which has been estimated in detail for low energy, integrated fluxes |I(|, does not necessarily apply 
to high-energy or differential fluxes. Therefore, also the /x/e lepton event ratio might suffer of uncertainties larger 
than ±5% in some energy-angle bins. Our empirical estimate of such errors in the statistical analysis is detailed in 
Appendix B. More precise estimates of the relative v M and v e flux uncertainties are in progress |^,Q. A new, ab 
initio, fully three-dimensional calculation of the atmospheric neutrino flux jifj| is also expected to shed light on these 
issues. 


x We will often use the term “neutrino” loosely, to indicate both y’s and u’s. Of course, we properly distinguish v from v in 
the input fluxes and in the calculations of cross sections and oscillation probabilities. 
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FIG. 1. Distributions of neutrino fluxes in terms of the neutrino zenith angle O, for three representative energies ( E „ = 1, 
10, and 100 GeV). Upper panels: flux. Middle panels: r/ e +V e flux. Lower panels: +V ll )/(v e +V e ) flux ratio. Solid 

lines and dots refer to the calculations in [|j] (HKKM’95) and 0 (AGLS’96), respectively. In this figure, geomagnetic effects 
are not reported, and the fluxes are symmetric under cos0 —> — cos0. Notice the increase of the v^/ve ratio with energy and 
for | cos0| —> 1. 


Concerning the overall uncertainty of the theoretical neutrino flux normalization, it is usually estimated to be 
20-30%. Most of the uncertainty is associated to the primary flux of cosmic rays hitting the atmosphere. It is 
important, however, to allow also for energy-angle variations of such normalization (e.g., the SK analysis Q| allows 
the spectral index to vary within ±0.05). Our approach to this problem is detailed in Appendix B. Here we want to 
emphasize that valuable information about the overall flux normalization can be obtained from more precise cosmic 
ray data from balloon experiments such as BESS (El]], CAPRICE |[I|], and MASS2 |52|j53| l. The BESS experiment 
has recently reported a relatively low flux of cosmic primaries [ p0|. which, as we will see, might represent a serious 
problem for the oscillation interpretation of the SK data (see also |l7| , [48f ). On the other hand, the MASS2 experiment 
can also measure the flux of primary protons and secondary muons at the same time p4| ], and might thus provide 
soon an important calibration of the theoretical flux calculations Therefore, it is reasonable to expect, in a 

few years, a reduction and a better understanding of the overall neutrino flux uncertainty, with obvious benefits for 
the interpretation of the atmospheric v anomaly. 
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In order to obtain measurable quantities (e.g., the lepton zenith angle distributions), one has to make a convolution 
of the neutrino fluxes with the differential cross sections and detection efficiencies (see, e.g., Hi®)- We consider 
five zenith angle distributions of leptons: sub-GeV muons and electrons, multi-GeV muons and electrons, and upward 
through-going muons. Concerning the calculation of the first four distributions, we have used the same technique 
used in [M for the old Kamiokande multi-GeV distributions. This approach makes use of the energy distributions of 
parent interacting neutrinos J5(J]. Concerning upward through-going muons, we improve the approach used in |l7j by 
including the zenith dependence of the SK muon energy threshold as given in ; we also use the same SK choice 

for the parton structure functions (GRV94 j57j, available in |38|) and muon energy losses in the rock [^9 201. For all 


distributions (SG, MG, and UP/r), we obtain a good agreement with the corresponding distributions simulated by 
the SK collaboration, as reported in Appendix A (to which we refer the reader for further details). 

In short, we can compute five distributions of SK lepton events as a function of the zenith angle 9 , namely, sub-GeV 
/u-like and e-like events (5+5 bins), multi-GeV /r-like and e-like events (5+5 bins), and upward through-going muons 
(10 bins), for a total of 30 observables. [] Few other analyses report explicit calculations of sub-GeV and multi-GeV 
zenith distributions (see, e.g., |^,^l|) or upward-going muon distributions (see, e.g., |i2|]) in agreement with the SK 
simulations. Other authors perform detailed calculations but use a reduced zenith information, as that embedded, 
e.g., in the up-down lepton rate asymmetry p3| (see, e.g., JmEH]). 


B. SK Data: Total and differential rates 


The experimental data used in our analysis are reported in Table I and Table II, together with the corresponding 
expectations as taken from the SK MonteCarlo simulations. The numerical values have been graphically reduced from 
the plots in ||,|il| and thus may be subject to slight inaccuracies. 

Table I reports the zenith angle distributions of sub-GeV and multi-GeV events collected in the SK fiducial mass 
(22.5 kton) during 535 live days, for a total exposure of 33 kTy. Fully and partially contained multi-GeV muons 
have been summed. Only single-ring events are considered. The distributions are binned in five intervals of equal 
width in cos0, from cos 9 = — 1 (upward going leptons) to cos 9 = +1 (downward going leptons). The total number of 
events in the full solid angle is also given. The quoted uncertainties for the data points are statistical. The statistical 
uncertainties associated to the SK MonteCarlo expectations originate from the finite simulated exposure (10 years 
live time [jl7| ) .p| 

Table II reports the differential and total flux of upward through-going muons as a function of the zenith angle. 
Data errors are statistical only. In this case, there is no statistical error for the SK theoretical estimates, which are 
derived from a direct calculation |2(| and not from a MonteCarlo simulation. 

It is useful to display the information in Tables I and II in graphical form. To this purpose, we take the central 
values of the theoretical expectations in Tables I and II as “units of measure” in each bin. In other words, all /r and 
e event rates (either observed, or calculated in the presence of oscillations) are normalized to their standard (i.e., 
unoscillated) expectations /tq and eoj^j The following notations distinguish the various lepton samples: 


SGe = sub-GeV electrons , 

SG/i = sub-GeV muons , 

MGe = multi-GeV electrons , 

MG [i = multi-GeV muons , 

UP/i = upward through-going muons , 
MC = theory (no oscillation) . 


2 Our earlier calculations of event spectra for pre-SK experiments can be found in E3-0Q-0- 
Systematic errors, not reported in Table I, are discussed in Appendix B. 

4 This representation was introduced in [[fi) to show those features of the atmospheric v anomaly which are hidden in the g/e 
ratio and emerge only when fi and e rates are separated. 


4 





TABLE I. Super-Kamiokande sub-GeV and multi-GeV e-like and /x-like atmospheric neutrino data, compared with the 
corresponding MonteCarlo simulations in different zenith angle (9) bins [jlj]. The numerical values are graphically reduced from 
the plots in [fl||T>]| . Units: Number of events. Experimental exposure: 33 kTy, corresponding to 22.5 kton fiducial mass x 535 
live days. Simulated exposure: 10 yr jlTj]. Errors are statistical only. The Super-Kamiokande MonteCarlo simulations refer to 
HKKM’95 neutrino fluxes §]. 


Event 

Bin 

cos 9 

Observed 

MonteCarlo 

Observed 

MonteCarlo 

sample 

No. 

range 

e-like events 

e-like events 

/x-like events 

/x-like events 

Sub-GeV 

1 

[-1.0, -0.6] 

287+ 16.9 

209 ± 5.5 

182 ± 13.5 

326 + 6.9 


2 

[-0.6, —0.2] 

231 ± 15.2 

206 ± 5.5 

225 ± 15.0 

316 + 6.8 


3 

[-0.2, +0.2] 

259 + 16.1 

220 ± 5.7 

228 + 15.1 

307 + 6.7 


4 

[+0.2, +0.6] 

227+ 15.1 

216 + 5.6 

264 ± 16.2 

308 ± 6.7 


5 

[+0.6, +1.0] 

227+ 15.1 

198 + 5.4 

259 + 16.1 

317 + 6.8 


total 

[-1.0, +1.0] 

1231 + 35.1 

1049 ± 12.4 

1158 ± 34.0 

1574 ± 15.2 

Multi-GeV 

1 

[—1.0, —0.6] 

50 + 7.1 

37 + 2.3 

64 + 8.0 

114 + 4.1 


2 

[-0.6, -0.2] 

56 ± 7.5 

51 + 2.7 

75 ± 8.7 

132 ± 4.4 

3 

[-0.2, +0.2] 

70 ± 8.4 

62 ± 3.0 

136 + 11.7 

173 ± 5.0 

4 

[+0.2, +0.6] 

74 ± 8.6 

52 + 2.8 

142 ± 11.9 

139 ± 4.5 

5 

[+0.6, +1.0] 

40 ± 6.3 

34 + 2.2 

114 + 10.7 

111 + 4.0 

total 

[-1.0, +1.0] 

290 ± 17.0 

236 ± 5.9 

531 ± 23.0 

669 ± 9.9 


TABLE II. Super-Kamiokande 535 day data on upward through-going 

muon fluxes, compared with the corresponding theo- 

retical calculations in different zenith angle (9) bins 

jlfiT. Units: 10' 13 cm 

_2 s _1 sr _1 . The calculated muon fluxes, taken from 

PP 

], refer to HKKM’95 neutrino fluxes |j5[, GRV’94 DIS structure functions j57|, and Lohmann muon energy 

losses in the 

rock [ 

5lJ. Errors are statistical only. The numerical values are graphically reduced from the plots in |l8[|l9|. 


Bin 

cos 9 


Observed 

Theoretical 

No. 

range 


/x flux 

/x flux 

1 

[-1.0, 

-0.9] 


1.03 + 0.18 

1.25 

2 

[-0.9, 

—0.8] 


1.16 + 0.18 

1.38 

3 

[-0.8, 

—0.7] 


0.90 + 0.17 

1.46 

4 

[-0.7, 

-0.6] 


1.62 + 0.22 

1.57 

5 

[-0.6, 

-0.5] 


1.31 + 0.18 

1.67 

6 

[-0.5, 

—0.4] 


1.57 + 0.20 

1.78 

7 

[-0.4, 

-0.3] 


1.59 + 0.21 

1.93 

8 

[-0.3, 

—0.2] 


2.20 + 0.25 

2.18 

9 

[-0.2, 

-0.1] 


2.73 + 0.28 

2.52 

10 

[-0.1, 

-0.0] 


3.42 + 0.31 

3.03 

total 

[-1.0, 

-0.0] 


1.75 + 0.07 

1.88 
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SuperKamiokande, total lepton rates (33.0 kTy, 535 day) 
M/ Mo m/ Mo m/ Mo 


0.6 1 1.4 1.8 

—I-1-1- 

UP/A ~ 

PC/A ‘ 

MC 


0.2 0.6 1 1.4 


UP/i - 
PC/A — 

MCX1.2 


0.6 1 1.4 1.8 


UP/A * 

PC/a ± 

MCX0.8 



(/a— e osc. favored) 


(yti —t osc. favored) 


(z/ osc. disfavored) 


Theory I MC = SK MonteCarlo with HKKM’95 z/ fluxes 

(±30% fi-e normaliz. error, ±5% /A/e ratio error) 

Data I SG = SubGeV events 

MG(FC) = MultiGeV events (fully contained) 

PC/a = Partially contained muons 
UP/a = Upward-going muons 
(±1a statistical error) 

FIG. 2. Super-Kamiokande data on total lepton events, compared with their theoretical predictions. The numbers of e-like 
and /a- like events are normalized to the central values of the corresponding MonteCarlo (MC) simulations eo and /ao jfj, obtained 
with HKKM’95 fluxes pi] (as reported in Table I). Error bars of experimental data are statistical; slanted error bars of MC 
represent ±30% systematics in the common /io,eo normalization and ±5% systematics in the /Ao/eo ratio. Left panel: Default 
MC (favors v^ oscillations). Middle panel: MCxl.2 (favors —* v T oscillations). Right panel: MCx0.8 (disfavors v 

oscillations). See the text for details. 


Figure 2 compares theory and data for the total lepton rates. The MG/a data sample is further divided into fully 
contained (FC) and partially contained (PC) events (notice that MGe events are all FC). The total number of lepton 
events are displayed in the plane (/a//a o,e/eo), so that the standard expectations (MC) correspond to the point (1,1) 
for each data sample. The UP and PC muons have no electron counterpart and are shown in the single variable 
n/Ho (upper strips). We attach a ±30% common uncertainty to the MC muon and electron rates (large slanted error 
bar), and allow for a ±5% uncertainty in the ///e ratio (small slanted error bar). The three subfigures 2(a), 2(b), and 
2(c) correspond to three choices for the theoretical predictions (MC): (a) Standard MC expectations; (b) MC rates 
multiplied by 1.2; and (c) MC rates multiplied by 0.8. 

Figure 2(a) clearly shows that, with respect to the standard MC, SK observes a deficit of muons (stronger for 
low-energy SG /a data and weaker for high-energy UP/a data) and, at the same time, an excess of electrons (both in the 
SG and MG samples). In the hypothesis of neutrino oscillations, these indications would favor —> v e transitions, 
with a mass square difference low enough to give some energy dependence to the muon deficit. As far as total SK 
rates are concerned, this is a perfectly viable scenario. 


6 












SuperKamiokande (535d) 


+ Data ± 1 a (stat.) 

zenith distributions 


_Theory ± 30% (no osc.) 



FIG. 3. Super-Kamiokande distributions of lepton events in terms of the lepton zenith angle 9 (cos 9 = —1, 0,1 correspond to 
upgoing, horizontal, downgoing leptons). From left to right: sub-GeV electrons (SGe) and muons (SG/x), multi-GeV electrons 
(MGe) and muons (MG/x), and upward through-going muons UP/x. In each bin, the observed rate R is divided by the expected 
rate Ro in the absence of oscillations, as taken from Tables I and II. Therefore, in this plot Theory=l with ±30% normalization 
error, and the deviations of Data (dots) from the flat theoretical expectations show the zenith anomaly at glance. 


Figure 2(b) shows how the previous picture changes when one allows for an overall increase of the MC expectations 
(say, ±20%). The relative excess of electrons disappears, while the muon deficit is enhanced. This situation is 
consistent with —* v T transitions, which leave the electron rate unaltered. The current SK data link rather strongly 
t'/x —*• v T oscillations to an overall increase of the MC expectations: indeed, a good fit to the SG and MG data samples 
requires a MC “renormalization” by a factor ~ 1.16 Q]. Although this factor is acceptable at present, it might not be 
so in the future, should the MC predictions become more constrained. In particular, if the recent BESS indications 
[ 501 for a relatively low flux of cosmic primaries were confirmed, then one should rather decrease the MC expectations 

i@. 

Figure 2(c) shows the effect of a MC decrease by 20%. In this case one would observe no deficit of muons (and 
even an excess of UP/x events) and a ~ 40% excess of electrons, which cannot be obtained in any known oscillation 
scenario. 

The discussion of Fig. 2 shows that: (i) Theory (no oscillation) and data disagree, even allowing for a MC renor¬ 
malization; (ii) The oscillation interpretation depends sensitively on the size of the renormalization factor; (iii) If this 
factor turns out to be < 1, the oscillation hypothesis is jeopardized; and (iv) It is thus of the utmost importance to 
calibrate and constrain the theoretical neutrino flux calculations through cosmic ray balloon experiments such 
as BESS JHiJ, and especially through simultaneous measurements of primary and secondary charged particles as in 
the forthcoming CAPRICE and MASS2 analyses Q|. All this information would be lost if the popular “/x/e double 
ratio” were used. Another piece of information that would be hidden by the double ratio is the fact that, in Fig. 2, 
the SG and MG data points appear to be very close to each other in the (/x//xo,e/eo) plane, while it was not so in 
Kamiokande ||j|. We are, however, unable to trace the source of such a difference (which is independent of the MC 
normalization) between SK and Kamiokande. 

Since the total lepton rate information is subject to the above ambiguities, one hopes to learn more from differential 
rates and, in particular, from the zenith distributions of electrons and muons. These distributions are shown in Fig. 3 
where, again, the rates R have been normalized in each bin to the central values of their expectations Rq (from Tables I 
and II). Therefore, the no oscillation case corresponds to “theory = 1” in this figure, with an overall normalization 
error that we set at ±30%. Deviations of the data samples (dots with error bars) from the standard (flat) distribution 
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are then immediately recognizable. The electron samples (SGe and MGe) do not show any significant deviation from 
a flat shape, with the possible exception of a slight excess of upward-going (cosf? —> —1) SG electrons. On the other 
hand, all the muon samples show a significant slope in the zenith distributions, especially for multi-GeV dataj^j Most 
of this work is devoted to understand how well two- and three-flavor oscillations of active neutrinos can explain these 
features of the SK angular distributions. 

We mention that additional SK measurements can potentially corroborate the neutrino oscillation hypothe¬ 
sis, namely, the stopping-to-passing ratio of upward-going muons [|l8|,^9|, neutral-current enriched event samples 
fjT|[ r i(i 691, and azimuth (east-west) distributions of atmospheric v events Jl^,|7^|. Such preliminary data will be 


considered in a future work. 

Finally, the SK data themselves could be used for a self-calibration of the overall neutrino flux normalization. 
In particular, one should isolate a sample of high-energy, down-going leptons with directions close to the vertical, 
so that the corresponding parent neutrinos would be characterized by a pathlength, say, L < 50 km and by an 
energy E v > 10-20 GeV. Then, for a neutrino mass difference Am 2 smaller than 10 -2 eV 2 jl|, the oscillating phase 
a A m 2 L/E v would also be small, and the selected sample could be effectively considered as unoscillated, thus 
providing a model-independent constraint on the absolute lepton rate and on the neutrino flux. The present SK 
statistics for strictly down-going, high-energy leptons, is not yet adequate to such a calibration. We will come back 
to this issue in the following. 


C. CHOOZ results 


The CHOOZ experiment |g4| searches for possible V e disappearance by means of a detector placed at L ~ 1 km 
from two nuclear reactors with a total thermal power of 8.5 GW. With an average value of L/E v ~ 300 krn/GeV, it is 
able to explore the v e — > V e oscillation channel down to ~ 10~ 3 eV 2 in the neutrino mass square difference, improving 
by about an order of magnitude previous reactor limits |39| ] . The sensitivity to neutrino mixing is at the level of a few 
percent, being mainly limited by systematic uncertainties in the absolute reactor neutrino flux. The ratio of observed 
to expected neutrino events is 0.98 ±0.04 ±0.04, thus placing strong bounds on the electron flavor disappearance [ p4| . 
The CHOOZ limits have been recently retouched (weakened) 1 25 as a result of the unified approach to confidence 
level limits proposed in 0 - 

The impact of CHOOZ for atmospheric oscillation searches and for their interplay with solar neutrino oscillations 
[0| has been widely recognized (see, e.g., 00 [Z3| 0 ) . Earlier studies of the interplay between reactor, atmospheric, 
and solar neutrino experiments can be found in 

Given their importance, we have performed our own reanalysis of the CHOOZ data in order to make a proper 
SK±CHOOZ combination. We use the v e + p —> e + ± n cross section as in our previous works Q, and convolute 
it with the reactor neutrino energy spectrum |24f | in order to obtain the positron rate. The expected rate is then 
compared with the data [^4],^5| through a y 2 analysis. We have checked that, in the case of two-family oscillations, 
we obtain with good accuracy the exclusion limits shown in ^5|. Our CHOOZ reanalysis will be explicitly presented 
in Section IV. The CHOOZ bound counts as one additional constraint (the observed positron rate); therefore, the 
global SK±CHOOZ analysis represents a fit to 30±1 observables. 


III. THREE-FLAVOR FRAMEWORK AND TWO-FLAVOR SUBCASES 


In this Section we set the convention and notation used in the oscillation analysis, 
among active neutrinos:^] 


We consider three-flavor mixing 


with 



U e 2 
Uy 2 
U t2 


U e3 \ 

U,3 
Ur3 ) 



(1) 


s It must be said that, in general, one cannot expect very strong zenith deviations in the SG data distribution, since the 
neutrino-lepton scattering angles are typically large at low energies (60°, on average) and therefore the flux of leptons is more 
diffuse in the solid angle. 

6 Oscillations into sterile neutrinos (see, e.g., |6l|[T7| fif} and references therein) are not considered in this paper. 
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mass(zq, zz 2 , ^ 3 ) = (mi,m 2 ,m 3 ) . 

It is sometimes useful to parametrize the mixing matrix U a i in terms of three mixing angles, u, <p, and ip: 


( 2 ) 


CcfrCu S(j) 

U ai = | SipSfiCoj C-xjjStjj + C'ljjCijj SipC(p 

~ Cijj S (f)Cijj + S^jSu C'fiSfpS^j S^C U C^Ccf) 


(3) 


where c = cos, s = sin, and we have neglected a possible CP violating phase that, in any case, would be unobservable 
in our framework. The mixing angles (oj,<p,ip) are also indicated as (^ 12 , 6 * 13 , 6 * 23 ) in the literature. 

While three-flavor oscillation probabilities are trivial to be computed in vacuum (i.e., in the “atmospheric part” of 
the neutrino trajectory), refined calculations are needed to account also for matter effects in the Earth. As in our 


previous works 36|44|], we solve numerically the neutrino evolution equations for any neutrino trajectory, taking into 


account the corresponding electron density profile in the Earth. Our computer programs are designed to compute the 
neutrino and antineutrino oscillation probabilities for any possible choice of values for the neutrino masses (mi, m 2 , m 3 ) 
and mixing angles (u>,<p,ip). 


A. Neutrino masses 

A complete exploration of the three-flavor neutrino parameter space would be exceedingly complicated. Therefore, 
data-driven approximations are often used to simplify the analysis J42j. As in our previous works ||2|,|3(|[39],[hi| , we 
use the following hypothesis about neutrino square mass differences: 

|m| — mi| = 5m 2 -C m 2 = |?n 3 — m\\ , (4) 

i.e., we assume that one of the square mass differences (5m 2 ) is much smaller than the other (m 2 ), which is the one 
probed by atmospheric neutrino experiments (and accelerator or reactor experiments as well |p9|). The small square 
mass difference is then presumably associated to solar neutrino oscillations |3^ j. 

Notice that the above approximation involves squared mass differences and not the absolute masses (which cannot 
be probed in oscillation searches). In particular, Eq. (|]) simply states that there is a “lone neutrino” z/ 3 , and a 
“neutrino doublet” (zq, zq), the doublet mass splitting being much smaller than the mass gap with the lone neutrino. 
However, Eq. (Q) can be fulfilled with either m 3 > mp 2 or m 3 < mi j2 . These two cases are not entirely equivalent 
when matter effects are taken into account, as shown in |^ 6 | . However, the difference is hardly recognizable in the 
current atmospheric v phenomenology [|36| . |65| . For simplicity, in this paper we refer only to the case m 3 > mp 2 , i.e., 
to a “lone” neutrino zq being the heaviest one. 

As far as m 2 > 1CV 4 eV 2 and 5m 2 <C 10 ” 4 eV 2 , atmospheric neutrino oscillations depend effectively only on 
m 2 . However, for larger values of 5m 2 the approximation ([|) begins to fail, and subleading, <5m 2 -driven oscillations 
can affect the atmospheric v phenomenology. We will briefly comment on subleading effects in Sec. VI C. Earlier 
discussions of such effects in solar and atmospheric neutrinos can be found in ]44|. 


B. Neutrino mixing 

Under the approximation (|J) one can show that CP violating effects are unobservable, and that the angle 10 can be 
rotated away in the analysis of atmospheric neutrinos (see |3C| | and references therein). In other words, atmospheric 
v experiments do not probe the mixing u> = 8 1 2 associated with the quasi-degenerate doublet (z^i,z^ 2 ), but only the 
flavor composition of the lone state zq, 


^3 — bf e3 u e A U [j,3 t'fj, A Ur3 W (b) 

— zz e A C(j)(sip A Cijj zq) . (6) 

The neutrino oscillation probabilities in vacuum assume then the simple form: 

-P vac (^a v a ) = 1 - 4C/ 2 3 (1 — U 2 3 ) S , (7) 

P vac K - v p ) = 4 U 2 a3 C/| 3 S (a + p) , (8) 


where 
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S = sin 1.27 


7 m 2 [eV 2 ] • L [km] 
E v [GeV] 


(9) 


neutrino and antineutrino probabilities being equal. However, in matter P(y) ^ P(y), and the probabilities must be 
calculated numerically for the Earth density profile. For a constant density, they can still be calculated analytically 


One can make contact with the familiar two-flavor oscillation scenarios in three limiting cases: 


4> = o => v 3 = Vfj + v T (pure zz M v T osc.) , (10) 

7r 

= - =>■ v 3 = S 0 v e + C 0 Vy, (pure <->• v e osc.) , (11) 

ip = 0 => v 3 = s^Vf, + c^Vt (pure v e <-> v T osc.) , (12) 

with the corresponding, further identifications (valid only in the 2v cases): 

sin 2 29= 4, (13) 

sin 2 2 0 etl = 4s|4 , (14) 

sin 2 2 9 eT = Aslcl . (15) 


Finally, we remind that for pure v tJ <-> v T oscillations the physics is symmetric under the replacement ip —> n/2 — ip, 
due to the absence of matter effects. Such effects instead break the (vacuum) symmetry <p —> 7 r /2 — <p for pure <-> v e 
oscillations, for which the cases <p < 7 r /4 and <p > 7 t /4 are distinguishable.R In generic three-flavor cases, no specific 
symmetry exists in the presence of matter and the mixing angles (p and ip must be taken in their full range [0, 7 t/ 2 ]. 
A full account of the symmetry properties of the oscillation probability in 2v and ?>v cases, both in vacuum and in 
matter, can be found in Appendix C of ||3G|]. 


C. The atmospheric v parameter space 

As previously said, under the hypothesis (f|) the parameter space of atmospheric v's is spanned by 
(m 2 , U e 3 , [/ m3 , U t 3 ). At any fixed value of m 2 , the unitarity condition 

Ul 3 + u% + u 2 3 = 1 (16) 

can be embedded in a triangle graph whose corners represent the flavor eigenstates, while a generic point 

inside the triangle represents the “lone” mass eigenstate v 3 . By identifying the heights projected from v 3 with the 
square matrix elements t/ 2 3 , U 2 3 , and f7 2 3 , the unitarity condition is automatically satisfied for a unit height 
triangle |3S[ |. 

Fig. 4 shows the triangle graph as charted by the coordinates U 2 3 (upper panel) or ( cp,ip ) (lower panel). When i/ 3 
(the mass eigenstate) coincides with one of the corners (the flavor eigenstates), the no oscillation case is recovered. 
The sides correspond to pure 2v oscillations. Inner points in the triangle represent genuine 2>v oscillations. 


IV. TWO-FLAVOR ANALYSIS 

In this Section we study first how the theoretical zenith distributions are distorted, in the presence of two flavor 
oscillations, with respect to the “flat” expectations of Fig. 3. This introductory study does not involve numerical fits 
to the data, and helps to understand which features of either <-> v T or <-► v e oscillations may be responsible for 
the observed SK zenith distributions. We then fit the most recent SK data (33 kTy) using a y 2 statistics, and discuss 
the results. 


7 Alternatively, one can fix p> < 7 r /4 and consider the cases m 2 > 0 or m 2 < 0 | 6 l| 
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FIG. 4. Triangle graph representing the three-flavor mixing space of v 3 . In the upper panel the parameter space is spanned 
by the matrix elements [7 2 3 , U 2 3 , and [7 2 3 . When such elements are identified with the heights projected by a generic point 
(xx 3 ), the unitarity condition ]U [7 2 3 = 1 is automatically satisfied (for a triangle of unit height). Vertices, sides, and inner 
area correspond to no oscillation, two-flavor oscillations, and three-flavor oscillations, respectively. In the lower panel, the same 
parameter space is charted through the mixing angles <p = #13 and ip = 623- 


A. Zenith distributions for <-» v T oscillations 

Figure 5 shows, in the same format as in Fig. 3, our calculations of the five zenith distributions of atmospheric 
neutrino events. In this and in the following figures, the upper left box contains comments on the scenario, while 
the upper right box displays the selected values of (m 2 /eV 2 , U 2 3 , U 2 3 , U 2 3 ). In Fig. 5 we consider, in particular, pure 
v fJ <-> v T oscillations {U 2 3 = 0) with maximal mixing (sin 2 26 = 4U 2 3 U 2 3 = 1). Of course, the electron distributions 

are not affected by v T transitions, while the muon event rates are suppressed, especially for zenith angles 

approaching the vertical (cos(9 = —1, upward leptons), corresponding to longer average neutrino pathlengths. The 
prediction for m 2 = 10~ 3 eV 2 (dashed line) is in reasonable agreement with all the muon data samples (SG, MG, and 
UP). For m 2 = 10~ 2 eV 2 the expected rates of SG/x and UP^x are significantly suppressed, and for m 2 = 1CU 1 eV 2 
one approaches the limit of energy-averaged 2v oscillations, with a flat suppression of ~ 50%, which does not appear 
in agreement with the data. On the other hand, decreasing m 2 down to 10 -4 eV 2 (thick, solid line), one has almost 
“unoscillated” distributions for the high energy samples MG/x and UP/x, since the phase m 2 L/i% is small. At lower 
energies (i.e., SG/x events), however, this phase can still be large enough to distort the zenith distribution. 

Notice that in Fig. 5 the theoretical electron distributions SGe and MGe are always below the data points. Using 
the overall ±30% normalization freedom, one can imagine to “rescale up” all the five theoretical distributions (by, say, 
15-20%) to match the electron data. This upward shift would also alter the muon distributions at the same time, but 
one can easily realize that such distributions would still be in reasonable agreement with the data for m 2 = 10 -3 and 
10 -2 eV 2 . Therefore, at sin 2 26 ^ = 1 one expects an allowed range of m 2 around m 2 = 10 -3 - 10 -2 eV 2 , independently 
of the details of the statistical analysis. Values of m 2 outside this range do not agree with the muon data. 
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As a final comment to Fig. 5, we notice that the theoretical rate for downgoing multi-GeV muons (rightmost bin 
of the MGp sample) is practically identical to the unoscillated case for m 2 < 1CP 2 eV 2 , as also observed at the end 
of Sec. II B. Therefore, the SK data can potentially self-calibrate the absolute muon rate normalization independently 
of oscillations , provided that the total experimental error in the last MG/i bin is reduced to a few percent. 

In Fig. 6 we take to 2 fixed (at 3 x 10 -3 eV 2 ) and vary the z^-z/ T mixing (sin 2 20 MT = 4?7 2 3 t/ 2 3 = 1,0.75,0.36). 
The suppression of the muon rates increases with increasing mixing; however, there is no dramatic difference between 
sin 2 2 6nr = 1 and 0.75 (solid and dashed lines, respectively). A value as low as sin 2 2 d^ T = 0.36 is not in agreement 
with the SG and MG muon distributions, although it is still allowed by the UP/z sample (which does not place 
strong bounds on the mixing). If all the distributions were renormalized (scaled up) to match the SGe and MGe 
samples, small mixing values would be even more disfavored. Therefore, we expect the mixing angle to be in the 
range sin 2 2 0^ T = 0.8-1, independently on the details of the statistical analysis. 

We emphasize that, although in principle e-like events are not affected by <-► v T oscillations, the SGe and MGe 
samples affect indirectly the estimate of the mass-mixing parameters, since they drive the fit to higher values of 
the neutrino fluxes. Further experimental constraints on the overall neutrino flux normalization might have then a 
significant impact on the current estimates of m 2 and sin 2 2 0^ T - 


B. Zenith distributions for v [L <-> v R oscillations 


Fig. 7 is analogous to Fig. 5, but for <-> v e oscillations with maximal mixing (t/ 2 3 = 0 and sin 2 0 e/i = 4t/ 2 3 t/ 2 3 = 1). 
Matter effects are included. The expected SG and MG rates of electrons coming from below appear to be enhanced, 
due to Uffs oscillating into zz e ’s. The slope of the zenith distribution is stronger for MGe than for SGe, because: (i) 
The v^/v e flux ratio increases with energy, as observed in Section II A; and (ii) the “angular smearing” due to the 
different lepton and neutrino directions is more effective for SG events. On the other hand, the suppression of the 
muon rates is not as effective as for the <-> v T case in Fig. 5. In fact, now there are some v e 's oscillating back into 
zz^’s. Moreover, matter effects tend to suppress large-amplitude oscillations (when the mixing is maximal in vacuum, 
it can only be smaller in matter). In general, one has a too strong increase of electrons and a too weak suppression 
of muons, although this pattern may be in part improved by rescaling down all the theoretical curves. In this case, 
the distributions at m 2 = 10~ 3 and ICC 2 eV 2 can get in marginal agreement with all the data, while m 2 = 10 -4 eV 2 
is in any case excluded. Notice, however, that to 2 = 10~ 2 eV 2 is not allowed by CHOOZ pit. 


^Oi/ e with sin 2 2^ eM =1 
Effect of varying m 2 
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FIG. 7. Distortions of the zenith distributions induced by variations of m 2 (eV 2 ), for pure ry, <-> i/ e oscillations with maximal 
(/•V,z/ e ) mixing. 
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Figure 8 shows the effect of varying the <->• v e mixing, for m 2 = 8 x 10~ 4 eV 2 —a value safely below the CHOOZ 
bounds. It can be seen that variations of the mixing do not help much in reaching an agreement with the data for 
this value of to , which thus seems to be disfavored. 

Figure 9 shows the size of matter effects for m 2 = 8 x 10~“ 4 eV 2 and maximal <-» v e mixing. For values of m 2 
around 10~ 3 eV 2 such as this, matter effects are more important for the multi-GeV sample, while lower m 2 values 
would enhance the effect for the sub-GeV sample and higher values for the UP/x sample (see, e.g., and Figs. 4 and 
6 in It can be seen that, in any case, the net matter effect is a decrease of the slope of the zenith distribution 

for both muons and electrons, i.e., a suppression of the oscillation amplitude with respect to the “pure vacuum” case. 
The effect is not completely reduced to zero above horizon, due to the neutrino-lepton angular smearing. 

We summarize the content of Figs. 7-9 by observing that the theoretical zenith distributions, in the presence of 
Vp <-> v e oscillations, are at most in marginal agreement with the SK data set. The agreement is somewhat improved 
by rescaling down the expectations. Earth matter effects are sizable and cannot be neglected in the analysis. 


C. Fits to the data 

In the previous two subsections we have presented quantitative calculations of the zenith distributions in selected 
two-flavor scanarios, and a qualitative comparison with the data. Here we discuss the results of a quantitative 2v fit 
to the SK data. 

Figure 10 shows the results of our x 2 analysis of the SK data (SGe, SG/x, MGe, MG/j, and UP/x data combined). 
The panel (a) refers to <-> v T oscillations (<f> = 0) in the plane (to 2 , sin 2 2ip). We find a minimum value Xmin = 29-6 
for 28 degrees of freedom (30 data points minus 2 oscillation parameters), indicating a good fit to the data.^| This 
is to be contrasted to the value Xmin = 126 for the no oscillation case, which is therefore excluded by the SK data 
with very high confidence. The quantitative limits on the mass-mixing parameters are consistent with the qualitative 
expectations discussed in Section IV A. Moreover, the allowed region is in agreement with the global analysis of 
pre-SK data shown in Fig. 2 of ||36| JF| Our allowed range of m 2 in Fig. 10(a) is somewhat narrower than the range 
estimated by the SK Collaboration jj. We have checked that the differences are largely due to the fact that only SG 
and MG were fitted in j|, while here we include also UP/x data, which help to exclude the lowest values of m 2 (see 
also |Q). To a lesser extent, our different definition of x 2 (see Appendix B) also plays a role. 

Figures 10(b) and 10(c) refer to <-> v e oscillations (ip = 7 t/2 ) in the plane (to 2 , sin 2 2 cp), for cp < 7 t/4 and (f> > 7r/4, 
respectively (the two cases being different, see Sec. Ill B). The minimum value of \ 2 is now much higher (67.7 and 68.6 
for panel (b) and (c), respectively) indicating that v M <-> v e oscillations are disfavored by the SK data. This represents 
an important step forward with repect to pre-SK data, which did not distinguish significantly between zx /; v r and 
Vf! <-> v e oscillations |Q (and, actually, showed a slight preference for the latter (|^J3^]). In Figs. 10(b,c), the C.L. 
limits around the minimum appear to be shifted to higher values of to 2 if compared to Fig. 10(a), as a consequence 
of matter effects that suppress the effective mixing in the lowest range of m 2 (see also p6|). 

The <-> v e allowed regions in Fig. 10(b) and 10(c) have a relatively scarce interest. On the one hand, they 
represent a poor fit to the SK data themselves. On the other hand, they are excluded by the CHOOZ reactor 
experiment. Figure bFd) shows the CHOOZ bounds, as derived by our own x 2 reanalysis, in good agreement with 
the limits shown in p5| . Such bounds are in contradiction with the v^ <-> v e allowed regions in Fig. 10(b,c) at 90% 
C.L., (although there might be a marginal agreement at 99% C.L.); therefore, we do not make any attempt to combine 
SK+CHOOZ data in a 2v analysis. 

Summarizing, our results for two-flavor oscillations are consistent with the 2v analysis of the SK Collaboration [Q , 
namely: (i) The no oscillation hypothesis is rejected with high confidence; and (ii) v T oscillations are largely 

preferred over <-> v e . On our part, we add the following nontrivial statement: (iii) The present SK bounds on the 
i/fj, <-> v T mass-mixing parameters are in good agreement with those obtained from the global analysis of pre-SK data 
in j3(| (including NUSEX, Frejus, IMB, and Kamiokande sub-GeV and multi-GeV data). 


8 The best-fit value of m 2 is not very meaningful at present. We prefer to focus on confidence intervals. 

9 In J3^] we obtained an allowed range for m 2 larger than the one reported by the Kamiokande Collaboration 0, presumably 
as a result of a different approach to the statistical analysis [ ^i||35|] . See also the comments of [^o| about the Kamiokande 
bounds in 0. 
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Two-flavor analysis 






FIG. 10. Two-flavor oscillation fits to the SK zenith distributions (SG, MG, and UP fi combined), (a) Fit for <-> v T 
(4> = 0) in the plane (m 2 ,sin 2 2 ip). The cases ip < n/4 and ip > 7t/4 are equivalent, (b) Fit for <-> v e (ip = n/2) in the plane 
(m 2 , sin 2 2(p), for (p < n/4. (c) Fit for <-> v e (ip = n/2) in the plane (m 2 , sin 2 20), for (p > 7r/4. The cases (b) and (c) are 
different, due to earth matter effects. The limits coming from the CHOOZ experiment are also shown in panel (d), as derived 
by our reanalysis. The solid and dotted curves correspond to 90% and 99% C.L., i.e., to variations of \ 2 — Xmin = 4.61, 9.21 
for two degrees of freedom (the oscillation parameters). 
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Finally, we comment upon recent claims j|8lf| of inconsistencies within the SK data, under the hypothesis of v tJ v T 
oscillations. The argument goes as follows. If <-> v T oscillations are assumed, the electron excess has to be adjusted 
by rescaling up the theoretical predictions. Then the overall muon deficit becomes close to maximal (~ 50%) for the 
SG and MG data (see our Fig. 2, middle panel). Such large suppression of the muon rate seems to suggest energy- 
averaged oscillations (i.e., large m 2 ). On the other hand, the zenith distortions require energy-dependent oscillations 
(i.e., relatively small m 2 ). Some difference between the “total rate” and the “shape” information, although exacerbated 
by the semiquantitative calculations of [pl| , indeed exists but is not entirely new, as it has already been investigated 
by the SK Collaboration [Q. In fact, the slide No. 18 of jis) shows the separate <-> v T fits to the “total rate” and 
“shape” data, the former preferring values of m 2 tipically higher than the latter. However, in the same slide |l^] one 
can see a reassuring, large overlap between the two allowed ranges (at 90% C.L.) for to 2 ~ few x 10~ 3 eV 2 . This means 
that, within the present (relatively large) experimental and theoretical uncertainties, there is no real contradiction 
between different pieces of SK data, as also confirmed by our good «-> v T global fit to the SK data. However, the 
above remarks, as well as our comments on the absolute event rates in Sec. II, should be kept in mind when new, 
more accurate experimental or theoretical information will become available. 


V. THREE-FLAVOR ANALYSIS 

In this Section we discuss in detail a three-flavor analysis of the SK and CHOOZ data. In the first subsection we 
show representative examples of 3v oscillation effects on the zenith distributions. In the second subsection we discuss 
some issues related to the L/E v variable. In the third and fourth subsections we report the results of detailed fits to 
the SK data, without and with the additional constraints from the CHOOZ experiment, respectively. 

We remind that, in three flavors, the CHOOZ mixing parameter sin 2 2 0 ee can be identified with 4Z7 2 3 (1 — U 2 3 ) 
[Eq. (j7|)]; therefore, the CHOOZ constraint sin 2 29 ee < 0.22, valid for m 2 > 2 x 10” 3 eV 2 p5f , translates into either 
U 2 3 < 0.06 or C / 2 3 > 0.94. 


A. Zenith distributions: Expectations for three-flavor oscillations 

In a three-flavor language, two-flavor oscillations with maximal «-> v T mixing are characterized by (U 2 3 , U 2 3 , U 2 3 = 
0,1/2,1/2) (the center of the lower side in the triangle graph of Fig. 4). Analogously, maximal <-> v e mixing is 
characterized by (U 2 3 , t/ 2 3 , U 2 3 = 1/2,1/2, 0) (the center of the right side of the triangle). A smooth 3v interpolation 
between these 2v cases can be performed by gradually increasing the value of U 2 3 from 0 to 1/2, and decreasing the 
value of U 2 3 from 1/2 to 0 at the same time, the element U 2 3 being adjusted to preserve unitarity. 

This exercise is performed in Fig. 11 for a relatively high value of m 2 (m 2 = 8 x 10 -2 eV 2 ). The thick and thin 
solid lines represent the zenith distributions for pure v T and pure <-> v e oscillations with maximal mixing, 

respectively. The dashed and dotted lines represent intermediate 3v cases, the first being “close” to v T (with an 

additional 20% admixture of v e ), and the second being “close” to v. L <-> v e (with an additional 20% admixture of v T ). 
Although none of the four cases depicted in Fig. 11 represents a good fit to all the SK data, and three of them are 
excluded by CHOOZ, much can be learned from a qualitative understanding of the zenith distributions in this figure. 

In the presence of oscillations, the distributions R/Rq in Fig. 11 are roughly given by 

— ~Pw + — Pe», (17) 

Mo Mo 

- ~ Pee + — P„e , (18) 

eo eo 

/ro and eo being the unoscillated rates. For a relatively high value of m 2 as that in Fig. 11, the asymptotic regime of 
energy-averaged oscillations approximately applies, except for the rightmost bins of the MG and UP/i distributions 
(where the shorter pathlengths require higher m 2 ’s for reaching such regime). Then, for pure <-» v T oscillations 
with maximal mixing, one has P ee = 1, P e ^ = 0, and ~ 1/2, so that e/eo = 1 and m/mo ~ 1/2, as indicated by 
the thick, solid lines in Fig. 11. 

For pure v e oscillations with maximal mixing, one has P ee ~ 1/2, P efl ~ 1/2, and P ~ 1/2, so that 

m/mo ~ (l + eo/Mo)/2 and e/eo ~ (l + Mo/eo)/2. For sub-GeV data, the often-quoted value Mo/eo ~ 2 applies, so that 
e/eo ~ 1.5 and \ij hq ~ 0.75, as indicated by the thin, solid lines in the SGe and SG [i panels of Fig. 11. For multi-GeV 
data, however, the value Mo/eo is not constant, ranging from ~ 2 along the horizontal (cos 9 = 0) to ~ 3 along the 
vertical (cos 6 = ±1, see also Fig. 1). Therefore, the ratio m/mo decreases slightly from ~ 0.75 (horizontal) to ~ 0.67 
(vertical), while the ratio e/eo increases significantly from ~ 1.5 (horizontal) to ~ 2 (vertical). This is particularly 
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evident as a “convexity” of the MGe distribution in Fig. 11 (thin, solid line). Such behavior is not confined to pure 
two-flavor oscillations, but is also present in genuine three-flavor cases, as indicated by the dashed and dotted lines 
in the MGe panel of Fig. 11. This shows that, in the presence of (v^,v e ) mixing (U 2 3 yf 0), the variations of the 
unoscillated ratio /xo/eo with cos 6 induce distortions of the zenith distributions even in the regime of energy-averaged 
oscillations |j^], contrary to naive expectations. Notice that such distortions do not depend on L/E Vl but on L and 
E v separately through no/eo = E v ). 

Another peculiar distortion, not dependent on L/E v , is related to a genuine three-flavor effect in matter This 

effect is basically due to the splitting of the quasi-degenerate doublet (zq, zq) in matter which, in the limit of large m 2 , 
leads to an effective square mass difference Sm 2 t oc E v and thus to a subleading oscillation phase Sni^^L / E u oc L 
which does not depend on E v but only on L |3q] . The main effect, relevant for atmospheric neutrinos, is to decrease 
the survival probability P ^ by an amount SP which, for a constant electron density N e , reads p8[ 


SP = 4 


U 2 e3 U 2 3 U: 


(i - U 2 3) 


r 3 • 2 
sin 


2.47(1 - U 2 3 ) 


N P 


rol/cm 


■ cos 0 


(19) 


Notice that the oscillation amplitude can be sizable only for large three-flavor mixing, while it disappears for two-flavor 
mixing (i.e., when one of the U 2 3 is zero), as indicated by the comparison of the dotted and thin solid lines in the 
MGfi and UP/x panels of Fig. 11. The phase of SP (the argument of sin 2 ) can be rather large in the Earth matter 
(N e ~ 2-6 mol/cm 3 ) and is modulated by the neutrino zenith angle 0. This modulation is particularly evident in 
the genuine 3zx cases of the UP/x panel in Fig. 11 (dotted and dashed lines), since upward through-going muons are 
highly correlated in direction with the parent neutrinos (9 ~ 0). The cos© modulation of SP is increasingly smeared 
out in the lower energy MG and SG muon samples. 
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FIG. 11. Three-flavor oscillations at rri = 8 x 10 


Distortions of the zenith distributions induced by variations of the 


2>v mixing. The values of U &3 interpolate smoothly between the pure 2 v subcases. Some cases are excluded by CHOOZ. 
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We remind that in Fig. 11 the dashed curves (U 2 3 = 0.2) correspond to a 3^ “perturbation” of pure <-> v T 
oscillations (thick solid curves, U 2 3 = 0), while the dotted curves (U 2 3 = 0.2) correspond to a 2>v “perturbation” of 
pure v^ <-> i/ e oscillations (thin solid curves, U 2 3 = 0). By comparing the dashed and thick solid curves, it can be 
seen that 3z/ perturbations of <-> v T oscillations do not alter dramatically the zenith distributions. The opposite 
happens for the v^ <-> v e case (dotted and thin solid curves). This pattern can be explained as an interference between 
vacuum and matter effects. 

More precisely, let us consider the case of large to 2 and no/s-o ~ 2. A3 v perturbation of <-» v T oscillations with 
maximal mixing can be parametrized by taking U 2 3 = e and U 2 3 = U 2 3 = (1 — e)/2 (the 2v case being recovered for 
e = 0). The relevant oscillation probabilities are then P ee ~ 1 — 2e, P ell ~ e, and P MAI ~ 1/2 — SP [notice that SP 
is of 0(e), see Eq. ©]• Equations (0® give the muon and electron rates, n/l^o ~ 1/2 + e/2 — SP and e/eo ~ 1, 
respectively. Therefore, to the first order in e = U 2 3 , the electron rate does not vary, and the muon rate varies little 
since e/2 and SP partly cancel. 

Conversely, a 3z^ perturbation of <-> v e oscillations with maximal mixing can be parametrized by taking U 2 3 = e 
and U 2 3 = U 2 3 = (1 — e)/2 (the 2v case being recovered for e = 0). The relevant oscillation probabilities are then 
P ee ~ 1/2, P ei j, ~ 1/2 — e, and P ~ 1/2 — SP. The muon and electron rates are now given by n/lM) ~ 3/4 — e/2 — SP 
and e/eo ~ 3/2 — 2e, respectively. To the first order in e = t/ 2 3 , the electron rate decreases, and also the muon rate 
is suppressed, since the terms e and SP have the same sign. Therefore, adding some v T mixing (U 2 3 ^ 0) to <-> v e 
oscillations changes the predictions considerably (generally in the direction of a better fit to the data). For this 
reason, we expect significant changes in the fit to SK data when moving continuously from pure <-> v e oscillations 
to genuine 3z/ cases. 

For simplicity, we have discussed the above three-flavor effects at large m 2 . Values of to 2 lower than in Fig. 11 are 
more interesting phenomenologically (being less constrained by CHOOZ) but more difficult to understand qualitatively, 
since the oscillations are no longer energy-averaged, and 2z/, 3z/, vacuum, and matter effects are entangled. Numerical 
calculations are required, and the results for m 2 = 8 x 1CP 3 and 8 x 10 -4 eV 2 are shown in Figs. 12 and 13, respectively 
(with mixing values chosen as in Fig. 11). The value m 2 = 8 x 10 -4 eV 2 is safely below the CHOOZ bounds. In 
Figs. 12 and 13, the genuine 3zz cases (dashed and dotted curves) show an improved agreement with the data (with 
respect to pure <-> v T oscillations), since they give an excess of electrons without perturbing too much the muon 
distributions. Although this advantage is not decisive at present, in view of the large uncertainties affecting the 
absolute normalization of the lepton rates, it might become crucial when such uncertainties will be reduced. 


B. Is L/E v a good variable? 

If two-flavor v tl <-> v T oscillations (U 2 3 = 0) were the true and exclusive explanation of the SK atmospheric data, 
then it would make sense to try to reconstruct the (unobservable) L/E v distribution of parent neutrinos from the 
lepton energies and directions. In fact, any distortion effect should be related to vacuum oscillations and thus to the 
ratio L/E v rather than to L and E u separately. The theoretical and experimental L/E v distributions for SK can be 
found in 0. 

However, beyond the 2v approximation, there are several oscillation effects that do not depend on L/E„. Some of 
these effects, originating from v e mixing ( U 2 3 ^ 0), have been described in the previous subsection. Non -L/E v effects 
also arise in the presence of two comparable square mass differences (i.e., 5m 2 ~ m 2 instead of 5m 2 <C m 2 ), or if 
non-oscillatory phenomena contribute to partially explain the data. 

Therefore, plots in the L/E v variable convey correct and unbiased information only under the hypothesis of pure 
Vp <-> v T oscillations. In all other cases, including our 3z^ framework, such plots cannot be used consistently. The 
following 3z/ fits, as for the 2v cases, make use of the zenith distributions and not of the reduced L/E v information. 
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C. Fit to Super-Kamiokande data 


Figure 14 shows our three-flavor fit to the SK data in the triangle graph, for values of m 2 decreasing from 2.5 x 10~ 2 
to 4.0 x 10 -3 eV 2 . The curves at 90% and 99% C.L. correspond to an increase of y 2 by 6.25 and 11.36 above the 
global minimum. The results, shown separately for sub-GeV, multi-GeV, and upward-going muons in the first three 
columns of triangles, are then combined in the last column. The CHOOZ data are excluded in this fit, in order to 
study what one can learn just from the SK data. 

For m 2 = 2.5 x 10 -2 eV 2 , the sub-GeV data exclude all two-flavor oscillation subcases (the triangle sides), but 
are consistent with genuine three-flavor oscillations at large U 2 3 (> 0.5). Also multi-GeV data are not in agreement 
with two-flavor oscillations, although the v T subcase (lower side of the triangle) cannot be excluded at 99% 

C.L. Multi-GeV data are well fitted by genuine three-flavor oscillations, but in a range of U 2 3 different (lower) than 
for sub-GeV data. The quality of the MG fits improves rapidly as one moves from the right side (z/ M <-> v e ) to the 
inner part of the triangle, as expected from the discussion of Fig. 11 in Sec. V A. The good 3v fit to SG and MG 
data is mainly driven by the genuinely 3zz matter effects discussed in Sec. V A [see Eq. (|l9|) and related comments]. 
Upward going muon data are much less constraining—at 99% C.L. they allow any oscillation scenario. At 90% C.L. 
they disfavor: (i) Pure or quasi-pure v e «-> v T and —> v e oscillations (left and right sides); and (ii) Large three- 
flavor mixing (the central region of the triangle). Large 3v mixing is excluded because it suppresses and distorts too 
much the UP/r distribution (see Fig. 11). The regions allowed separately by SG, MG, and UP data have no common 
intersection, and the combination of all the data is a null region (last triangle). 

For m 2 = 1.5 x 10~ 2 eV 2 , the increasing energy-angle dependence of the oscillation probability helps to fit the data 
better. Therefore, the regions allowed in the triangle are larger for each of the three data sets, and a small allowed 
region appears at 99% C.L. in the combination. Such region is enlarged for a lower value of m 2 (1.0 x 10 -2 eV 2 ) and it 
appears also at 90% C.L. at to 2 = 6.5 x 10~ 3 eV 2 . For the latter value of to 2 , all the data are consistent with <-> v T 
oscillations (the lower side), while the other two-flavor subcases (the left and right sides) are excluded. The MG 
data sample “repels” v t , <-> v e oscillations more strongly than SG data since, as observed in Sec. II A, higher-energy 
samples are characterized by a larger /i/e unoscillated ratio, and thus are more sensitive to the presence (or absence) 
of mixing. 

For to 2 ~ 4.0 x 10~ 3 eV 2 , also the global combination of the data is consistent with v IL <-> v T oscillations, with large 
but not necessarily maximal mixing. Large values of U 2 3 are also allowed, indicating that the SK data, by themselves, 
do not exclude large three-flavor mixing. We remark that the goodness of the fit improves rapidly when one moves 
from the right side inwards (i.e, when one “perturbs” pure <-> v e oscillations), while it changes more slowly when 
one moves from the lower side upwards (i.e., when one “perturbs” pure v {1 <-> v T oscillations), as expected from the 
discussion in Sec. V A. 

Figure 15 is analogous to Fig. 14, but for lower values of m 2 , ranging from 2.5 x 1CV 3 to 4.0 x 10 -4 eV 2 . As to 2 
decreases, the SG data fit is not affected very much, since values as low as fewxlO" 4 eV 2 still provide a good fit to 
this sample. However, the oscillation phase starts decreasing more rapidly for higher-energy samples (MG and UP), 
leading to an insufficient suppression of the muon rates and to a gradual reduction of the allowed regions. Notice, in 
particular, how UP/i data constrain 3zz mixing for to 2 ~ 1-3 eV 2 . In any case, the preferred regions are more and 
more reduced and closer to the lower side of the triangle, corresponding to smaller allowed values for U 2 3 . There is 
no joint allowed region at 99% C.L. for iri 2 below ~5x 10~ 4 eV 2 . 

We summarize the three-flavor fit to SK data as follows: (i) The SK data exclude both U 2 3 = 0 and U 2 3 = 0, being 
consistent with large <-> v T mixing (not necessarily maximal); (ii) Values of U 2 3 as large as 0.5 cannot be excluded 
only on the basis of SK data. Indeed, 3zz oscillations with large u e mixing can improve the fit and, in particular, it can 
explain (part of) the electron excess in the SG and MG samples. Quantitative bounds on the mixing matrix elements 
I/ 2 3 can be derived from Figs. 14 and 15. In the next subsection we study the impact of CHOOZ on such indications. 
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FIG. 14. Three-flavor analysis in the triangle plot. Separate and combined fits to the Super-Kamiokande atmospheric 
neutrino data, for selected values of m 2 . The parameter space is defined in Fig. 4. Notice that the allowed regions for the 
combined fit are always close to the lower side of the triangle, i.e., to pure <-> v T oscillations, although it is not necessarily 
so for the separate data samples. Relatively large values of mixing ( U 2 3 ) are allowed. 
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D. Fit to Super-Kamiokande and CHOOZ data 


We combine SK and CHOOZ data (30+1 observables) through a joint % 2 analysis. The results are reported in 
Fig. 16. 

Figure 16 shows the 90% and 99% C.L. in the triangle plot for selected values of m 2 , ranging from 4 x 10 -3 eV 2 
(upper triangles) to 6.5 x 10 -4 eV 2 (lower triangles). The left column of triangles reports the fit to SK data only (as 
derived from Figs. 14 and 15). The middle column reports the fit to CHOOZ data, which exclude a large horizontal 
stripe. In fact, the nonobservation of v e disappearance implies that v 3 is either very close to the upper corner v e 
(so as to suppress oscillations) or very close to the lower side {v^ <-► v T oscillations being unobservables in CHOOZ). 
Clearly, the addition of the CHOOZ bounds to the SK fit (right column of triangles) cuts significantly the upper 
part of the solutions, so that only relatively low values of U 2 3 are allowed. However, the CHOOZ bounds are rapidly 
weakened as m 2 is decreased, and for m 2 = 1.5 x 10 -3 eV 2 the parameter U 2 3 can be as large as ~ 0.15 at 90% 
C.L. and ~ 0.25 at 99% C.L., corresponding to a significant v e appearance probability [see Eq. (||)]. Therefore, the 
CHOOZ data constrain but do not exclude the role of v e mixing and electron appearance in the interpretation of the 
SK data (see also |||,H). In particular, part of the electron excess in the SG and MG samples could be explained 
by nonzero values of UJ 3 rather than by uncertainties in the overall neutrino flux normalization. Nonzero values of 
U 2 3 also contribute to distort the zenith distributions |l8j, as discussed in Sec. V A. 

So far we have seen the impact of CHOOZ on the mixing parameters U 2 3 . The impact on the square mass difference 
to 2 is summarized in Fig. 17, which shows the x 2 as a function of to 2 , for unconstrained values of the mixing angles 
and for both fits to SK (dashed lines) and SK+CHOOZ (solid lines). The minimum value of y 2 is 28.3 for SK and 
29.8 for SK+CHOOZ, indicating a good fit to the data (30 and 31 observables, respectively). The CHOOZ data help 
to constrain to 2 on the higher range, but its role decreases rapidly for to 2 < 10 -3 eV 2 . 

Also shown in Fig. 17 are the 90% and 99% C.L. intervals for to 2 , which allow values as low as 5 x 10~ 4 eV 2 . 
The possibility of exploring such low values of to 2 should be seriously considered in long baseline experiments. An 
interesting result of Fig. 17 is the stability of the to 2 range indicated by SK—it does not change dramatically by 
adding the CHOOZ constraint. Therefore, the inclusion of the CHOOZ data in the global analysis affects more the 
mixing than the mass parameter. 

A comparison of SK and pre-SK bounds is illuminating. Figure 17 should be compared with Fig. 10 of [p6| , where we 
combined the data from NUSEX, Frejus, IMB, and Kamiokande (sub-GeV and multi-GeV) in a three-flavor analysis. 
The comparison shows that the SK and SK+CHOOZ bounds on to 2 are perfectly consistent with the pre-SK bounds. 
The SK+CHOOZ data appear to improve significantly the old upper bound on to 2 , but give a lower bound very similar 
to the pre-SK data. Notice that we have long since claimed that the popular value to 2 ~ 10~ 2 eV 2 overestimated the 
best fit for pre-SK data, and that values as low as 5 x 10 -4 eV 2 were compatible with the atmospheric v data |^6| . 
We plan to perform a joint analysis of all the data (SK+CHOOZ+pre-SK) in a future work. 


VI. IMPLICATIONS OF THE 2,v ANALYSIS 

In this Section we examine some implications of our three-flavor analysis for the phenomenology of atmospheric, 
long-baseline, and solar neutrino experiments, as well as for model building. 


A. Atmospheric v phenomenology 

The SK atmospheric data are consistent with 3z^ oscillations with dominant v^ <-> v T transitions (U 2 3 ■ U 2 3 > 0.2 at 
90% C.L.) and subdominant v e <-> z/^. T mixing ( U 2 3 < 0.15). The mass square difference to 2 is favored in the range 
~ 0.8 x (10~ 3 -10~ 2 ) eV 2 . Can one improve significantly such indications only with SK or other atmospheric data? 

Large v >J; v T mixing should generate a v T flux comparable to the flux at the detector site. However, the 
“contamination” of /u-likc and e-like events from r production and subsequent leptonic decay is estimated to be very 
small in SK Q. Therefore, there seems to be little hope to test the v T channel through r appearance in SK. 

Nevertheless, the possibility of enhancing (through appropriate cuts) the r —> /r and r —> e “pollution” in selected 
and e event samples may deserve further attention in other atmospheric v detectors such as Soudan2. 

The tests of v e mixing (i.e., of the matrix element U 2 3 ) can certainly be improved with higher statistics SK data, 
in particular with more multi-GeV upgoing electron and muon events. Such data samples are characterized by a 
relatively high z+/z+ flux ratio (see Sec. II A), and thus are more sensitive to an increase of e-like events due to 
<-> z+ transitions. Multi-GeV data are already more powerful than SG (and UP/r) data in constraining U 2 3 (see 
Fig. 15). 
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FIG. 16. Three-flavor analysis in the triangle plot. Separate and combined fits to Super-Kamiokande and CHOOZ. The 
CHOOZ data exclude large horizontal stripes in the triangle plots. The combined SK+CHOOZ solutions are closer to pure 
i'll <-> v T oscillations as compared with the fit to SK data only. However, the allowed values of U 2 3 are never negligible, especially 
in the lower range of m 2 . 
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FIG. 17. Bounds on m 2 for unconstrained mixing, as derived from the y 2 analysis of SK data, with and without CHOOZ. 
The 90% and 99% C.L. intervals correspond to variations of \ 2 — Xmin = 6.25, 11.36 for three degrees of freedom (the oscillation 
parameters). 


While the elements t/ 2 3 determine the amplitude of oscillations, which can be already derived from total event rates, 
the parameter nr governs the phase of oscillations, and thus it can be derived only through event spectra. Hypothetical 
spectra of neutrino events as a function of E v and L would be the most sensitive probes of to 2 . Unfortunately, a 
complete kinematical closure of {/-induced events cannot be achieved in SK, so neither E v nor L can be precisely 
reconstructed, especially for low-energy events. This intrinsic feature will eventually limit the maximum accuracy of 
to 2 fits attainable with SK data only. In this respect, the possibility of improving the E v and L reconstruction in 
experiments as Soudan2 |Q (through observation of the struck nucleon), or in high-density detectors as proposed in 
pdj, appears extremely interesting and promising. 
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B. Long Baseline experiments 


Long baseline (LBL) accelerator experiments, such as K2K j25|, MINOS |28|, and various CERN-Gran Sasso 
proposals [^0| , are expected to confirm the atmospheric neutrino signal with a controlled beam. Since both two-flavor 
[Q and three-flavor analyses like ours show that to 2 can be as low as ~ 5 x 10 -4 (99% C.L.), the design of low-energy 
beams should be pursued seriously. If the atmospheric m 2 range can be covered completely, then it suffices to have 
either a disappearance or a v T appearance signal to confirm the SK anomaly. 

However, we think that long baseline experiments should be designed to measure oscillation parameters, rather 
than merely to confirm an oscillation effect already found by SK. Measuring the oscillation parameters is a task that 
demands careful considerations, especially if 2>v oscillations are to be tested (for 2v oscillations, see the lucid discussion 
in H]). 

The determination of m 2 requires energy spectra analyses, and thus high-statistics event samples. If m 2 happens 
to be in the low range of the experimental sensitivity, the r appearance sample might consist of just a handful of 
events. The e appearance event sample might also be small if Uf 3 —> 0. Therefore, a safe reconstruction of w 2 should 
be based mainly on fj, event spectra from the v tJ —> disappearance channel, where most of the signal is expected in 
any case. This implies a good monitoring of the initial is fl beam with a near detector. 

The determination of the matrix elements U 2 3 requires that several oscillation channels are probed at the same 
time—redundancy is never enough to constrain neutrino mixing p0[ . For instance, the disappearance channel 

is sensitive only to U 2 3 but tells nothing on Uf 3 or Uf 3 , while the z/ ;i —> v T appearance channel is sensitive to the 
product U 2 3 U 2 3 but it cannot separate the two factors U 2 3 and Uf 3 nor measure U 2 3 [see Eqs. (j7j.||)]. These aspects 
of ?>v mixing tests ||^,^] in long-baseline experiments are better appreciated in Fig. 18. 

Figure 18 shows, in the first column of triangles, the region allowed by SK+CHOOZ for selected values of to 2 . The 
second column shows, superimposed, the prospective regions that can be probed at ~ 90% C.L. by K2K Q, both in 
the v fl —> channel (slanted bands) and in the z+ —> v e channel (hyperbola).^] It appears that K2K might not reach, 
in the v tl —► v e channel, sufficient sensitivity to probe the values of U 2 3 allowed by SK+CHOOZ. The disappearance 
channel can cover the whole SK+CHOOZ region, but only for m 2 > 2 x 10~ 3 eV 2 . Therefore, K2K is 

basically expected to give information on U 2 3 for m 2 >2x 10 -3 eV 2 , given the sensitivities prospected in p9| . 

With respect to K2K, the MINOS experiment is being designed to probe lower values of m 2 and to explore also the 
-*■ v t appearance channel (the region below the hyperbola touching the —> v T side in Fig. 18). Possible signals 

in the three channels —► v e ,n,T will constrain the quantities U 2 3 U 2 3 , U 2 3 , and U 2 3 Uf 3 , respectively, so that the 

elements U 2 3 can be pinpointed for to 2 > 2 x 10~ 3 eV 2 if the uncertainties are kept small. For lower to 2 ’s, MINOS 
rapidly looses sensitivity in at least one of the oscillation channels, and it might be difficult to constrain the neutrino 
mixing parameters. 

Notice that, for U 2 3 ~ 0.15 and the U 2 3 ~ 0.5 (allowed by SK+CHOOZ), the /+ —► v e appearance probability is 
Pfj,e = ^U 2 3 U 2 3 (S) ~ 0.3 • (5), where (S)(< 1) is the oscillation factor in Eq. (|), averaged over the v beam energy 
spectrum. Depending on (S) and on the specific mixing parameters, values of P Me as large as 15% appear possible in 
properly designed LBL experiments. 

A final remark is in order. The sensitivity regions in Fig. 18 have been derived from the prospective estimates 
reported in the experiment proposals [pf)[[28[ , which are in continuous evolution (even more so for the CERN to Gran 
Sasso proposals [Q, not shown). Therefore, the above considerations on K2K and MINOS are to be considered 
as preliminary and qualitative. Nevertheless, it remains true that LBL experiments might face some difficulties in 
constraining the mixing parameters, especially if m 2 is low or if the three oscillation channels cannot all 

be probed. Re-directing the goal of LBL experiments from “confirming the Super-Kamiokande signal” to “measuring 
the parameters (m 2 , Uf 3 , U 2 3 , Uf 3 )” would be beneficial to the current debate on various LBL proposals. 


10 Curves of isoprobability are either of the form (7 2 3 = const in the v a —> v a disappearance channel or of the form U 2 3 Up 3 = 
const in the v a —> vp appearance channel 
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FIG. 19. Beyond the approximation of one mass scale dominance: example of effects induced by the subdominant mass scale 
5m 2 , for oscillation parameters allowed by atmospheric, CHOOZ, and solar neutrino data. 


C. Solar neutrino problem 

In the limit 5m 2 -C m 2 [Eq. (^|)], experiments with terrestrial (atmospheric, accelerator and reactor) neutrino 
beams probe the parameters (to 2 , U 2 3 , U 2 3 , U 2 3 ). In the same approximation, solar neutrino experiments probe the 
parameters (5m 2 , U 2 X , U 2 2 , U 2 3 ), i.e., the small mass square difference and the mass composition of u e [ fi2| . 

Therefore, both terrestrial and solar v experiments probe the element U 2 3 , which parametrizes the mixing of i/ e with 
the “lone” neutrino mass eigenstate u 3 . The effect of such mixing on solar neutrinos is to give an energy-independent 
contribution to the disappearance of zz e ’s from the sun. Since this contribution is basically proportional to the square 
of U 2 3 H^], sizable mixing is required to have large effects. However, as U 2 3 is constrained by the SK+CHOOZ fit, 
such effects are relatively small to be detected in the current solar neutrino experiments. The smallness of U 2 3 also 
reduces the “coupling” of the terrestrial and solar parameter spaces Q. 

Nevertheless, it is interesting to notice that a preference for small values of U 2 3 emerges naturally from solar neutrino 
data only, in both matter-enhanced [^3|,^2) and vacuum |86) three-flavor oscillation fits. An updated analysis using the 
latest SK solar neutrino data would be desirable to confirm such indication. In fact, a more pronounced preference of 
solar i/’s for relatively small values of U 2 3 would be a nontrivial, although indirect, hint that both solar and terrestrial 
v data are consistent within the same oscillation scenario. 

The relation between solar and atmospheric neutrinos is not necessarily confined to a preference for relatively small 
values of the mixing matrix element U 2 3 . For instance, if the assumption in Eq. (||) is violated, atmospheric neutrinos 
can become sensitive to the subleading oscillations driven by 5m 2 , i.e., by the “solar neutrino” mass difference. 
Figure 19 shows a representative case beyond the 5m 2 <SC to 2 approximation. The values chosen for the oscillation 
parameters (5m 2 , m 2 ,ui, are allowed at 99% C.L. by the present SK+CHOOZ bounds and by the pre-SK, three- 
flavor solar v analysis in 132]. The effect of the subdominant mass scale 5m 2 can be evaluated by comparing the 
solid curves (5m 2 = 0) and the dashed curves (5m 2 yf 0). Of course, the effect is more significant when 5m 2 L/E v is 
larger, i.e., for upgoing SG events, while it rapidly decreases at higher energies (MG and UP/z samples) and shorter 
pathlengths (cos 6 —» 1). For SG/z events, however, the effect approaches the size of the error bars and thus it might 
be probed by SK! However, the strategies for disentangling the oscillations driven by to 2 and 5m 2 are nontrivial and 
will be investigated in a future work. 
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D. Models of neutrino mass and mixing 


Theoretical or phenomenological models of neutrino mass and mixing try to predict, or “explain,” the set of 
parameters (mi, m 2 , m 3 ) and U a i. Our analysis constrains the subset of parameters (m 2 ,U 2 3 ), provided that Sm 2 
is sufficiently small (< 1CP 4 ). Many models that try to explain solar+atmospheric v data fall within this category, 
and are thus strongly constrained by the SK+CHOOZ bounds worked out in this paper. For instance, the so- 
called bimaximal mixing model (?(J, characterized by (U 2 3 ,U 2 3 ,U 2 3 ) = (0,1/2,1/2) for atmospheric neutrinos [and 
by (t/ 2 ^ t/ 2 2 , U/ 3 ) = (1/2, 1/2,0) for solar neutrinos] is allowed for m 2 > 10 -3 eV 2 (see Fig. 15). Conversely, the 
trimaximal mixing model |87|], characterized by (U 2 3 ,U 2 3 ,U 2 3 ) = (1/3,1/3,1/3) for atmospheric neutrinos [and by 
(U 2 3 , U 2 2 , U 2 3 ) = (1/3,1/3,1/3) for solar neutrinos] appears to be strongly disfavored by our combined SK+CHOOZ 
analysis (it would correspond to the center of each triangle in Figs. 14 and 15). Of course, many other models (see, 
e.g., the classification in J88|) can be tested through our bounds on the oscillation parameters, provided that the 
dominance of the mass scale m 2 is assumed for atmospheric neutrinos. 

Other models try to explain the atmospheric anomaly and the LSND evidence | p9| for —► v e oscillations, with 
an allowance for an energy-averaged suppression of the solar neutrino flux. Arguments disfavoring such scenarios 
are discussed in Q (see, in particular, Table VI of Q and related comments). One such model has been recently 
proposed in [^o|, where to 2 ~ 0.4 eV 2 is assumed to drive the oscillations in the LSND experiment |8tJ range, as well 
as energy-averaged oscillations of atmospheric Vs, while Sm 2 ~ 10 -4 10~ 3 eV 2 is assumed to drive energy-dependent 
oscillations of atmospheric Vs. Both Sm 2 and m 2 can then contribute to the solar neutrino deficit through energy- 
averaged oscillations. Since the bounds worked out in Sec. V D assume Sm 2 <C 10” 4 eV 2 , and thus do not apply to 
such model, we have performed a numerical analysis of SK data ad hoc, using the same mass-mixing parameters as 
in pc)] . We find that the resulting zenith angle distributions of muons (not shown) are only mildly distorted, and 
that the model is disfavored by the SK atmospheric data at > 99% C.L., with or without matter effects. We mention 
that, for the choice of parameters in |po|| , matter effects influence significantly the zenith distributions, making them 
flatter than in vacuum. The semi-quantitative calculations in showed a more optimistic agreement to the SK 
data, in part because matter effects were ignored. In addition, the large value for U 2 3 chosen in J9G] does not appear 
in agreement with the global analysis of laboratory neutrino oscillation searches (including the LSND data) performed 
in joij]. 


VII. SUMMARY AND CONCLUSIONS 

We have performed a three-flavor analysis of the SK atmospheric neutrino data, in a framework characterized by 
the mass-mixing parameters (m 2 ,U 2 3 ,U 2 3 ,U 2 3 ), in the hypothesis of one mass scale dominance. The variations of 
the zenith distributions of v events in the presence of flavor oscillations have been investigated in detail. Fits to the 
SK data, with and without the additional CHOOZ data, strongly constrain the parameter space. Detailed bounds 
have been shown in triangle graphs, embedding the unitarity condition U 2 3 + U 2 3 + U 2 3 = 1. The allowed regions 
include the subcase U 2 3 = 0, corresponding to pure +, <-> v T oscillations. However, values of U 2 3 > 0 are also allowed. 
In particular, for m 2 close to (or slightly below) 10~ 3 eV 2 , U 2 3 can be as large as ~ 0.15 (at 90% C.L.). Scenarios 
with U 2 3 > 0 correspond to genuine three-flavor oscillations and are characterized by a rich phenomenology, not 
only for atmospheric Vs, but also for solar and laboratory neutrino oscillation searches. In particular, challenging 
opportunities are disclosed for v e appearance searches in long baseline experiments. Our analysis also places strong 
constraints on models of neutrino mass and mixing. In addition, we have examined many facets of the SK data and 
of their interpretation, that will deserve further attention when the experimental and theoretical uncertainties will be 
reduced. 
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APPENDIX A: CALCULATION OF ZENITH DISTRIBUTIONS 


The calculation of the zenith angle distributions of SG and MG lepton events involves the numerical evaluation of 
multiple integrals of the form 


$ <gi a ® e ® P , 


(Al) 


where $ is the unoscillated neutrino spectrum, er is the differential cross section for lepton production, e is the detector 
efficiency for lepton reconstruction, and P is the oscillation probability (h||3(J. 

The efficiency function e is not always reported in the experimental papers. In particular, it has not been explicitly 
given by the SK Collaboration so far. We faced a similar problem in the analysis of the Kamiokande multi-GeV data 
performed in pd| |. Our solution [ f36| was to use the energy distribution of the parent neutrinos v that produce a 
detected lepton £, namely, 


dN v 

dE u 


d$ v 

dE v 


J dE 4 


(A2) 


where Eg, is the lepton energy. This distribution, which gives information on the factor •!> <8> a ® e in Eq. (|aT|) , has 
been published in Q for the Kamiokande experiment. Concerning SK, we have used the analogous information from 


Using the energy distribution of parent neutrinos, it is possible to reconstruct the zenith distribution of the final 
leptons, provided that the smearing induced by neutrino-lepton scattering angle is taken into account p^ |. While 
for the old Kamiokande multi-GeV data we approximated this effect with an energy-independent smearing angle of 
~ 17°, for SK we properly take into account the distribution of the lepton scattering angle and its dependence on the 
energy, which is especially relevant for SG events (92|. We find good agreement with the SK estimate of the average 
scattering angle as a function of energy (as reported in |2l|], p. 99). Concerning the neutrino fluxes, we refer to || 
except for SG events, where we use the differential spectra from || with geomagnetic corrections p3| . 

Since the distributions of parent neutrinos in |p6| are given in arbitrary units, we need to normalize the total 
area of our estimated SK lepton distributions (SG and MG, in the absence of oscillations) to the corresponding 
values simulated by the SK Collaboration, as reported in Table I (total rates). For SG events, this renormalization 
compensates, in part, for the fact that we use low-energy v fluxes from HQ instead than from 0. Of course, a 
more direct calculation of the SG and MG distributions (avoiding the use of indirect information such as the parent 
v distributions) is preferable; we intend to perform such calculation when the SK efficiency function e(Eg) will be 
made publicly available. In any case, our present approach produces results in satisfactory agreement with SK zenith 
distributions for SG and MG events [Q, as shown in Fig. 20. The small differences between our calculations and the 
SK simulations are not relevant, being comparable to the SK MonteCarlo statistical error. 

Figure 20 also shows the UP/r distribution, for which we use a direct computation as in [ pTf, with the following 
ingredients: GRV94 DIS structure functions J57J, Lohmann et al. muon energy losses in the rock | |59|j20| , and the zenith 
dependence of the SK muon energy threshold from Also for this distribution, we obtain a good agreement 

with the corresponding SK calculation (with the same inputs). 

Notice that Fig. 20 refers to the no oscillation case. Some “oscillated” /i-like and e-like event distributions (as well 
as their ratio fi/e) have also been presented by the SK Collaboration in various Conferences, especially for the case 
of maximal <-> v T mixing. We obtain good agreement with SK also in such cases (not shown). 

In conclusion, we are confident that our calculations of the zenith distributions represent a satisfactory approxi¬ 
mation (not a substitute, of course) of the SK simulations. Improvements of our calculations for the SG amd MG 
samples are possible (with a more accurate knowledge of the SK detector efficiency) but do not appear to be decisive 
at present, in view of the good agreement reported in Fig. 20 and of the relatively large theoretical uncertainties 
discussed in the following Appendix. 
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(no oscillation) 



COS COS ^ COS 1? cos $ cos $ 

FIG. 20. Angular distribution shape: comparison of the SK MonteCarlo simulations (solid lines) with our calculations 
(dashed lines). 


APPENDIX B: STATISTICAL ANALYSIS 

Fitting histograms is a delicate task. In general, the predictions in any two bins are correlated, and ignoring such 
correlations typically leads to significant variations in the allowed ranges for the fit parameters. Concerning the 
SK zenith distributions, this problem adds to the difficulty of evaluating the size of uncertainties associated to the 
theoretical neutrino flux calculations and to the neutrino cross section, as well as the variations of such uncertainties 
in terms of the neutrino energy and direction. 

There is no easy solution to such problems, other than continually improving the calculations, understanding the 
role and the uncertainties of any input parameter, and removing as many approximations as possible [49|. Cosmic 
ray experiments can also help to constrain the models of atmospheric showers. The confidence in the estimated 
cross sections might benefit from a resurgence of experimental interest for low-energy neutrino interactions. In the 
meantime, it is wise to adopt conservative error estimates. 

For each of the 30 SK zenith bins {& ! :} j :=i,..., 3 o used in the analysis, we take the experimental statistical error Si from 
Tables I and II. The 30 x 30 statistical error matrix s?- simply reads s?- = SijSiSj. As a global systematic error a for 
each bin, we assume conservatively ±30% of the theoretical prediction (with or without oscillations). The systematic 
error matrix is then u?- = PyOjOj, where the correlation matrix py is evaluated as follows. 

If the systematic uncertainties Oj had a single, common origin such as an overall normalization uncertainty, then the 
bin values b-f's would be fully correlated (py = 1) and the systematics would cancel in any bin ratio bi/by However, 
the presence of several sources of uncertainties implies that py < 1 and that the ratio 6 ? ;/ b ;j is affected by a residual 
uncertainty 


a 2 (bi/bj) = of + of - 2 p ij a i a j , 

where all o’s represent fractional errors. For m = ctj = o(= 0.3) the above relation can be inverted to give 


Pij — 1 


2 (bi/bj) 

2 o 2 


(Bl) 


(B2) 
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which allows to estimate p^ from the ratio error. For instance, if bi refers to downgoing SG e-like events and 
bj to the corresponding xx-like events, with a u/e uncertainty of, say, ±5%, the corresponding correlation index is 
Pij ■ = 1 - (0.05 1 2 3 4 5 6 7 8 9 10 11 )/(2 ■ 0.30 2 ) = 0.986 || . 

The task is then reduced to the evaluation of the most important sources of errors for the ratios bi/bj. The 
total error for the /z/e flavor ratio (including the theoretical uncertainties and the experimental misidentification) is 
conservatively estimated to be ± 8 % for SG events and ±12% for MG events in [Q. For bins of equal flavor, one 
expects an additional energy-dependent uncertainty in the ratio bi/bj due to uncertainties in the neutrino energy 
spectrum slope. In fact, by comparing the relative rates of SG, MG, and UP/z events calculated with different input 
fluxes (either J(| or ||), we find typical ratio errors of ±5% for 6j(SG)/6j(MG) and 6j(MG)/6j(UP), and of ±10% 
for &i(SG)/6j(UP), i.e., errors increasing with the relative difference between the mean energies of the event samples. 
Finally, one expects also angular-dependent errors for b, 
between | cos 6 i | and | cos 6 

Qualitatively, all this means that the correlation between any two bins decreases from unity as the bins are more 
separated in energy, angle, and flavor, thus giving to the theoretical distributions some freedom to vary their shape. 
Quantitatively, we formalize the above estimates by generalizing Eq. (B2) as 


/ bj , that we estimate to be at most ±5% when the difference 


j,j is maximal (ratio of vertical to horizontal direction bins). 


Pij — 1 


2<r 2 


u e 

2 CT 2 


2a 2 ’ 


(B3) 


where er = 30%, and: (i) 07 is the “flavor-dependent uncertainty,” equal to 10% for bins of different flavors and zero 
otherwise; (ii) <je is the “energy-dependent uncertainty,” equal to zero for bins (i,j) belonging to the same sample 
(SG, MG, or UP), to 5% for bins (i,j) of the kind (SG,MG) or (MG,UP), and to 10% for bins (i,j) of the kind 
(SG,UP); and (iii) erg is the “direction-dependent uncertainty,” equal to 5% times the difference between the mean 
direction cosines |(cos0)i| and | (cos 0)y |. For instance, the first bin of the SGe distribution and the last bin of the UP/x 
distribution have the lowest correlation, ~ 0.874, since they are the most distant in energy, flavor, and direction. 

We finally define our y 2 function as 


\ 2 - Abi {slj + pijCTiCTj) 1 A bj , 
ij 


(B4) 


where A b is the difference between the bin contents in Tables I and II and our theoretical calculations (with or without 
oscillations). We mention that the y 2 fit to the SK data appears to be rather sensitive to 07 . Lowering its value from 
our present choice (10%, comparable to the estimates in Table II of [Q) to a few percent would shrink significantly the 
allowed regions but would also worsen the best fit. A reduction of this and other systematics would greatly improve 
the statistical power of oscillation hypothesis tests. 
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